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Motivation 


Source  localization  map  (SLM) 


Range  * 


*  Depth 


Localization  of  acoustic  sources  using 
passive  sonar 

Fundamental  task  for  monitoring  and 
surveillance  systems 

Difficult  due  to  complexities  of  the 
propagation  environment 

•  Ocean  behaves  as  an  acoustic  waveguide 

•  Varying  sound-speed  causes  acoustic  signals 
to  bend 

*  Temperature,  pressure  and  salinity 
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Modeling  Preliminaries 


Measurements 
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Our  data: 

{y/M  e  cw,v/}‘_i 


Acoustic  model 
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Model-based 

predictor 
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UW  acoustic  model 


Pfc,/(rfc) e  CJ 


N 


Replica 


Problem  statement 

Given  {y/M  e  c*v,v/H=1  an  UW  acoustic 
model  and  k,  estimate  the  source 
locations  {r*(i)}f=i  while  being  robust  to 
model  uncertainty. _ 
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Related  Works 


Matched-field  processing  (MFP)  [Baggeroer  ’88],  [Csenzak  ’97],  [Debever’  07],  [Mantzel  ‘12],  [Forero  ‘13] 


•  Measurements  y/(f)and  UW  acoustic  model 

•  Grid  of  tentative  locations  G  ■■=  {r9}°=1 

•  Use  UW  acoustic  model  to  generate  replicas  P/  ==  |pi(/(ri),...,pG>/(r0)l 

•  “Match”  measurements  and  replicas  to  construct  ambiguity  surface  - 


Acoustic  power 
estimates 


Matched-field  tracking  (MFT)  [Bucker  '94],  [Wilmut  ’98],  [Fialkowsky  '01] 


Score  all  possible  tracks 
across  ambiguity  surfaces 

Unconstrained:  Ot 
Constrained:  G  •  BT~^ 


Caveats:  . 


Computationally  taxing  due  to  batch  processing 

Temporal  information  is  not  used  when  constructing  ambiguity  surfaces 

Multiple  sources 


A  Nonlinear  Model  for  the  Measurements 


Data  model: 


K 

y f(t)  =  ^2sk,f(t)  [pkj(rk)+vkj(t,rk)]  +  ef(t),  V/ 

fc=i  | 

Acoustic  gain 

Fourier  coefficient  Perturbation 

for  jfcrth  source  at  ^  f 


Goal:  find  estimates  for  {rfc(t)}^_isequentially 

Caveats:  -  Mapping  h  ■.  rfc  is  not  available  in  closed  form 

-Perturbations  {vfcJ  foreland  acoustic  gains  (sfcj(f)}  are  unknown 
-  How  to  summarize  the  history (y/(r)  e  cw,v/}‘“\  ? 
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Challenges  Faced  by  Kalman-Based  Trackers 


Assumption:  v*j(^rk)  =0,  Vt3rfe 


Measurement  equation: 


Problematic! 


K 


y/(*)  =  Xlsfe,;(0p*j(r*)  +  e/(*)>  v/ 


a=i 


State  equation: 


x(f)  =  A x(t  -  1)  +  C{*) 


Acoustic  gains,  source  location, 
kinematics 


Sparsity-based  methods:  [Maliutov  ‘05],  [Edelman  ‘11],  [Fannjiang  ‘10],  [Mantzel  ‘12],  [Liu  ’12],  [Forero  '13, ’15] 


Grid:  Q  :=  {r9}°=1 

.  Depth 


◄ - 

Range 


New  model:  a 

y f(t)  =  ^sgj(t)pgj  +  e/(t),  V/ 


Large  dimensionality 
can  be  problematic 


s/(t)  =  As/(t  -  1)  +  C(t),  V/ 

I 

Sparse 


Sparse  Kalman  trackers  [Filos  ’13],  [Farahmand  ’14]  - ►  Performs  covariance  matrix  updates 
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Coherent  Broadband  Source  Localization 


Broadband  SLM 


Depth 


G 


New  model:  y f(t)  =  +  e/(i),  V/  S 

<j— 1 


1  Range 

|i  A  ||  ^ 

r«i  IK9II2 


GxF 


Coherence:  Common  support  across  ^  Q  ‘s 


S(t)  =  argmin  t  j^||y/(i)-P/s/||^  +  x  ^2  l|s/  -  s,(t-  1)|||  +  ||«9||2 

sgfcxF  i  t-:  &  ..  — 


/=! 


/= 1 


9= 1 


Summary  term  for  (y/(r)  e  C jV,  V/}*  =\ 


Remark:  Perturbations  can  be  captured  in  this  estimator  as  in  [Forero  ‘15] 
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Equivalent  Problem  and  Proximal  Gradient 


Real-valued  problem: 


-f1  "Lt  Y.T 

S(t)  =  ajgmin  )-  ^  ||y/(t)  -  P/S/|||  4-  |  ||vff  -  v3(f  -  1)||1  +  ||v„||2 

^  j=l 


=1  I 

V,  :=  Z'g+aY  €  R2F 


Proximal  Gradient  (PG)  method:  A  form  of  majorization-minimization 

Differentiable  part  of  the  cost  — __________  G 


min  A(S)  +  fi  ||  v. 


<J=1 


Majorizer:  function  H  is  a  majorizer  for  h  iff  #  (si  Z)  >  /i(S) 

H(  Z;  Z)  =  h{  Z) 

Lipschitz  constant  for  Vh(S) 

F  l 

PG  uses:  ff(S;Z):=fc(z)+^vfc/(«/),(a/-*/)+-^||S-z||5. 

/=! 
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PG  Solvers 


For  each  t,  iteratively  solve 


SW  (t)  —  argiiiin 

L. 

s 


G 


o-l 


Equivalent  to 


S^^argmiii^My  wg  -w],t_11(i)  2+^||v< 
S  ' 


Closed  for  solutions  per  (<) 


^  Decomposable 
across  g 


1  - 


+ 


Convergence:  For  any  p  >  o,  the  sequence  {S  converges  to  S (t)  as  i  -¥  oo, 

featuring  worst-case  convergence  rate  0{\/i). 

Proof:  Follows  from  Beck  ‘09] 


SWell EX-3  Dataset:  MFP 
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•  10  frequencies  between  53  Hz  and  197  Hz 

•  N  =  9  and  G  =  20, 000  (100  depths  and  200  ranges) 

•  M  =  lfor  SLM  and  M  =  10  for  all  others 

•  Per  time  frame,  plot  100  largest  peaks 


Single  frequency  149  Hz 
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Depth  (m)  ~  Depth  (m) 


SWellEX-3  Dataset:  Broadband  SLM  Tracks 


Constraint  region: 

22  depths  by  5  ranges 

Scores  averaged 
over  all  10  frequencies 

High  computational 
Complexity 

Roughly  77  min 
per  tracklet  (8  min) 
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Few  seconds 
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SWellEX-3  Dataset:  Multiple  Sources 
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Data  set  was 
artificially  modified 
to  emulate  two  sources 


MFT:  Challenging  to  identify  two 
different  tracks  due  to  width  of  the 
beams  in  the  ambiguity  surfaces 


Tuning  parameters:  grid  search  over  ^ 

=  mjx||[|pJilyi(i)  +  /isSli(i  -  1)1,  ■ .  ■ ,  |p*iFyf(*)  +  pSg .F(t  -  1)[ 
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Can  identify  the 
two  sources 


Selection  of  tuning 
parameter  A 
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Depth  (m) 


SWellEX-3  Dataset:  SLMs  and  Algorithm 


Two-source  case  at 
initial  position  (time  =  0) 


APG  solver  converges  much 
faster  than  PG  solver 


SLM  obtained  after  20 
APG  iterations 
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needed  to  identify  the  support  of 
the  SLM 
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Summary 


•  Sparsity  enables  underwater  source  tracking  of  multiple  sources 

•  Possible  to  coherently  process  broadband  data  for  improved  tracking 

•  Fast  PG  and  accelerated  proximal  gradient  (APG)  algorithms  for  constructing 
SLMs  and  source  tracks 

•  Reduced  computational  complexity  via  screening  of  predictors 

•  Tuning  parameter  selection  is  challenging 


Current  work: 

•  Validation  of  results  with  other  datasets 

•  Rules  for  selection  of  tuning  parameters 
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